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ABSTRACT 

We present a non-parametric method for decomposition of the light of disk 
galaxies into disk, bulge and bar components. We have developed and tested 
the method on a sample of 68 disk galaxies for which we have acquired I-band 
photometry. The separation of disk and bar light relies on the single assumption 
that the bar is a straight feature with a different ellipticity and position angle from 
that of the projected disk. We here present the basic method, but recognise that 
it can be significantly refined. We identify bars in only 47% of the more nearly 
face-on galaxies in our sample. The fraction of light in the bar has a broad range 
from 1.3% to 40% of the total galaxy light. If low- luminosity galaxies have more 
dominant halos, and if halos contribute to bar stability, the luminosity functions 
of barred and unbarred galaxies should differ markedly; while our sample is small, 
we find only a slight difference of low significance. 

Subject headings: galaxies: fundamental parameters — galaxies: photometry 
- galaxies: structure — galaxies: spiral — galaxies: kinematics & dynamics - 
methods: data analysis 



1. Introduction 



Hubble (1926) divided disk galaxies into barred and unbarred families, denoted SB and 
S respectively. A similar visual classification scheme was described by de Vaucouleurs (1959), 
who extended it to include an intermediate (SAB) family to allow for weak bars, although he 
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noted that the distribution of bar strengths seemed to be continuous. Sandage & Tamman 
(1981) designate intermediate types S/SB. Subjective classifications into two or three bins 
might be a useful first step, but a more quantitative estimator of the prominence of the bar 
would clearly be superior. 

Visual classification rests on the fact that the eye readily distinguishes at least three 
separate light components of a barred galaxy: the disk, bulge, and bar. Ideally, we would 
like a reliable procedure to separate the bar light from that of the disk and bulge in high- 
quality digital images of a large sample of galaxies. Such information would be useful to 
challenge models of the origin of bars in galaxies. More generally, the distributions of bar 
light fractions, shapes, and sizes in galaxies ought to be predicted by theories of galaxy 
formation and evolution. 

We will assume that the separation of the light into disk, bar, and bulge components 
is well founded. One worries whether the separate components are dynamically distinct; 
e.g. it is far from clear that every star pursues an orbit that is always confined to the one 
component. Nevertheless, it seems reasonable to hope that the amount of light in all three 
components does not change by much on an orbital time scale and that evolution is slow. 

We wish to measure a number of photometric properties of bars: the luminosity of the 
bar relative to that of the disk or of the entire galaxy, the length of the bar, the major axis 
light profile of the bar, the axis ratio or ellipticity of the bar, the shape of the bar, etc. Bars 
range from bright inner ovals (not very different from lenses) to skinny, more "boxy" features 
(Athanassoula et al. 1990), and from so-called "flat" bars, with a major-axis light profile that 
declines slowly, to "exponential" bars, that have a steeply declining light profile (Elmegreen 
& Elmegreen 1985). Such a multi-dimensional parameter space requires something more 
sophisticated than a single number to quantify the properties of a bar. 

Many attempts to characterize bars in galaxies have been reported; early work was 
summarized by Sellwood & Wilkinson (1993). Abraham et al. (1999) fit ellipses to the 
image at many isophote levels, and evaluate (&/ a )bar from deprojection of the inner ellipses 
assuming the outer ellipses indicate the projection geometry of a round, flat disk. Buta & 
Block (2001) define bar strength from the gradients of the gravitational potential in the 
disk plane deduced from a photometric image, rectified to face-on, and assuming a fixed M/L 
ratio and disk thickness. Laurikainen, Salo & Rautianen (2002) show that Qb correlates well 
with the maximum ellipticity of the bar. Jogee et al. (2004) and Marinova & Jogee (2006) 
also base their method on multiple ellipse fits, but define a galaxy to be barred from a 
combination of the ellipse parameters. Laurikainen, Salo & Buta (2005) quantify the non- 
axisymmetric light by Fourier decomposition. Radically different approaches were adopted 
by Seigar & James (1998), who quantified the bar flux after subtraction of an axisymmetric 
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disk, and by Prieto et al. (2001), Peng et al. (2002), and de Souza, Gadotti & dos Anjos 
(2004) who fit assumed parametric forms for the disk, bulge and bar to the photometric 
image. 

Abraham & Merrifield (2000) compare their bar-strength parameter with optical classi- 
fications of the Frei et al. (1996) sample of bright galaxies, finding some relation only when 
it is combined with a second parameter characterizing the concentration of the total light. 
Buta et al. (2005) estimate Qb for the Eskridge et al. (2000) sample of galaxies and find a 
fairly smooth distribution of values, suggesting a continuum of bar strengths from zero to 
a maximum. Marinova & Jogee (2006) reanalyse the same sample by their method, finding 
instead a preponderance of strong, i.e. quite elliptical, bars. 

Most of these studies aim for a single estimator, such as bar strength, and do not 
attempt to characterize the full multi-dimensional zoo of bars. In addition, many methods 
suffer from a number of intrinsic weaknesses. 

First, all require an estimate of the inclination and position angle of the projected 
disk; estimating projection geometry from the outer isophotes, assuming a flat, circular disk, 
can be complicated by the presence of non-axisymmetric features, such as the bar itself as 
well as spiral arms, oval distortions or warps in the outer disk. Recognizing where such 
features might be important and adjusting estimates appropriately requires practice and 
good judgement, processes that are intrinsically difficult to automate. 

Second, since the ellipticity at intermediate radii is affected by spiral arms and non- 
circular "rings" , as well as by bars, deciding the radius at which a bar ends and a different 
bisymmetric feature begins can be fraught with difficulty Recent discussions of the ap- 
propriate way to estimate bar lengths have been given by Athanassoula & Misiriotis (2001), 
Aguerri et al. (2003), Erwin (2005), Gadotti & de Souza (2006), etc. Buta, Block & Knapen 
(2003) describe a Fourier method to separate the bar from the spiral before computing Q b . 

Third, the ellipticity in the inner disk is affected by the bulge, which is generally a weak 
contributor to the light by the end of the bar, but usually has the highest surface brightness 
in the center. Furthermore, the Qb parameter requires that the bulge be subtracted from 
the disk and bar before deprojection and then added back (Laurikainen et al. 2004). 

Here we propose a new method to decompose the light of a galaxy in a high-quality 
image. Our method offers a more automated, and possibly superior, approach to all three 
of the above difficulties, and yields estimates of the bar light fraction, luminosity profile and 
ellipticity. We present a preliminary description here, but recognize that there is substantial 
scope for improvement. For example, we assume an elliptical light distribution for the bar, 
which is generally a poor fit; however, our method could be extended to include more general 
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shapes. 

Our technique is a generalization of the disk-bulge decomposition method described by 
Barnes & Sellwood (2003, hereafter Paper I). We assume that the disk is thin, flat, and round; 
the bar is also thin and flat, and is concentric with the disk. We base our decomposition of 
the bar and disk upon the assumption that the bar is straight and, in general, has a different 
position angle and ellipticity from that of the projected disk. We make no assumptions 
about the light profiles either of the disk or of the bar. Following Palunas & Williams 
(2001), Barnes & Sellwood also used a non-parametric form for the bulge; here we adopt 
instead a more conventional Sersic bulge model (Sersic 1968), but a parametric bulge model 
is not a required feature of our three-component decompositions. 

We reduce the subjectivity of defining the projection geometry by fitting the entire 
light distribution at once. Since a single elliptical isophote can be affected by local non- 
axisymmetric features, most workers generally look for a radial range over which the elliptic- 
ity and position angle do not vary much. By requiring a single ellipticity for the entire disk 
we evaluate a weighted average automatically. The average ellipticity and position angle 
in a disk-bulge decomposition can be thrown off by prominent non-axisymmetric features. 
Barnes & Sellwood attempted to estimate the systematic uncertainties caused by such fea- 
tures, by determining the spread in the fitted projection angles after the image is rectified 
and reprojected about other major axes. By fitting the most prominent non-axisymmetric 
feature, the bar, we significantly reduce this potential systematic error. The influence of 
bright spiral arms remains, but they should have a lesser effect both because they are gen- 
erally less strong than bars and because they, by their nature, do not have a fixed position 
angle. 

We avoid judging where the bar turns to a spiral and/or ring simply by fitting a straight 
feature of arbitrary light profile; the bar intensity profile should simply drop to the noise in 
the outer disk. Generally we find that the fit assigns some light to the bar beyond the end of 
the bar, as judged by eye, which is probably due in part at least to other non-axisymmetric 
features in the disk that are not completely orthogonal to the bar direction. As will be 
described below, we cut off the bar light profile at some distance from the center to prevent 
it from rising again where it crosses an outer spiral arm, for example. 

Finally, by fitting a separate bulge component, we are able to follow the disk and bar 
light profiles in towards the center, where they are required to have the same ellipticity and 
position angle as at larger radii. 

The data we use here were acquired specifically for this study. Our objective was to 
obtain high-quality images of a manageable, but representative, sample of galaxies selected 
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without regard to morphological appearance. We have used this sample to develop our 
technique and present the statistics of disk, bulge, and bar light fractions. We hope that the 
method can be refined and at least partly automated in order to be able to analyse the light 
of much larger galaxy samples, e.g. SDSS (York et al. 2000), the Millennium Galaxy survey 
(Liske et al. 2006), etc. 

2. Data 
2.1. Sample 

We have selected a magnitude-limited sample of galaxies with known redshifts in as 
unbiased a manner as possible. The redshift is important only to determine whether the 
properties of the fitted bar vary with absolute magnitude. Our starting point was the 
Southern Sky Redshift Survey (da Costa et al. 1998), which gives redshifts of a complete 
sample of galaxies selected from the Hubble Space Telescope Guide Star Catalog (Lasker et 
al. 1990) to a magnitude limit of m GSC < 15.5. Da Costa et al. estimate the uncertainty 
in the quoted magnitudes, which were determined automatically from scans of sky survey 
plates, to be about 0.31 magnitudes. 

We targeted a complete sub-sample of 129 objects from the SSRS in the declination 
range —16 < 5 < —14 and right ascension 10 hr < a < 16 hr. There were no large clusters 
of galaxies among the targets, but several small groups were apparent. 

2.2. Observations 

The galaxies were observed on the nights of April 6-10, 2000 using the CTIO 0.9 m 
telescope and the Tek2K#3 CCD. The CCD was operated unbinned, providing an image scale 
of 0.396"/pixel. Quad readout mode was used, with an approximate gain of 3 electron/ADU 
and read noise of 5 electrons per readout pixel. The standard CTIO B, V and I filters were 
used. Observing conditions were photometric throughout the first three nights, while the 
later part of the fourth night and the entire fifth night were plagued by varying amounts of 
cloud. Seeing conditions were variable throughout the run, with image FWHM varying from 
1" to 3", with 1.5" typical. Standard stars from the fields of Landolt (1992) were observed 
on the first four nights. Twilight sky flats were taken on the first, third, and fourth nights, 
and bias frames were taken at the beginning of each night. 

The observing scheme was designed to provide a deep, high-quality I band image to be 
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used to model the galaxy's light distribution. We also acquired less-deep images in the B 
and V bands in order to measure the colors of the modeled components. For each galaxy 
we took two B exposures, two V exposures, and three I exposures; the telescope was jogged 
slightly between exposures to minimize the effects of CCD artifacts. Exposure times for each 
image were 110 seconds in B, 75 seconds in V, and 240 seconds in I; on the last night the 
exposure times were doubled. 

We had time to observe 95 of the identified list of 129 target galaxies. The selection 
of galaxies to observe was based purely on their position, without reference to any intrinsic 
property. Most of the 34 omitted galaxies were in the more southerly half of the declination 
range. 

2.3. Reductions 

The images were overscan and bias subtracted, trimmed, and flattened with the twilight 
sky flats in the usual way, using IRAF0 We determined the centroids of stars in the images 
which we used to align multiple exposures; integer pixel translations were adequate, since 
the pixel scale was small compared to the seeing. We combined the aligned images with 
cosmic ray rejection. We estimated the sky brightness in the vicinity of each galaxy using 
the biweight (Beers et al. 1990) and subtracted it from the images. We also examined each 
image and masked out foreground stars and any remaining artifacts. 

We measured the brightnesses of standard stars in a fixed digital aperture of diameter 
2.5 times the worst seeing FWHM, and subtracted the local sky from a surrounding an- 
nulus. We did not attempt to measure the atmospheric extinction, but used the standard 
CTIO extinction coefficients from Stone Sz Baldwin (1983). We determined the photometric 
zeropoints for the first three nights of the run, with an internal precision of 4 to 5%. 

3. Modeling 

We attempt to divide the light of each galaxy into distinct components: a disk, a bar 
and a bulge. Our technique is an extension of that used for fitting a 2-D photometric image 
with a non-parametric disk and bulge model, which is described in the appendix of Paper 



1 IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the As- 
sociation of Universities for Research in Astronomy, Inc., under a cooperative agreement with the National 
Science foundation. 
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I. While Peng et al. (2002) and de Souza et al. (2004) also fit 2-D images, they assume the 
light can be decomposed into a number of components of fixed parametric form. Kent & 
Glaudell (1989) attempted a non-parametric fit to the light distribution of NGC 936, but 
modeled the bar along its major axis only. 

Here we make the usual assumption that the disk is thin, flat, and intrinsically round; 
the fitted ellipticity and position angle of the disk defines the projection geometry, therefore. 
The bulge is assumed to be spheroidal with its symmetry plane coincident with that of the 
disk, and we also require the disk, bulge and bar to be concentric. As in Paper I, we depart 
from convention by allowing the radial profile of the disk light to be non-parametric, but 
here parameterize the bulge as a spheroidal Sersic model. 

The new aspect of this work is that we allow for the possibility of a bar component in 
the disk, which we define to be an additional feature in the 2-D light distribution that is bi- 
symmetric with no spirality. Thus the disk and bar components are each constrained to have 
fixed ellipticity and position angle at all radii, but the light profiles of each have arbitrary 
tabulated values. As will be clear from our fits, our assumption of an elliptical shape for the 
bar light is generally poor; a more boxy shape would be more appropriate. Since we suspect 
that a more sophisticated characterization of the bar light (e.g. Athanassoula et al. 1990) 
would not have a strong effect on our conclusions, we leave this refinement for future work. 



Here {A} are the sky-subtracted intensities of the N pixels used in the image, with their 
associated uncertainties cr*. The photometric model is described by a tabulated set of inten- 
sities {Ik}, which are interpolated to the pixel locations using the weights w^i, and I s ^ is 
the contribution of the model bulge to the intensity of the ith pixel. 

We first describe the simplest case in some detail, which is to fit only a thin, circular 
disk that is inclined to the line of sight. The model isophotes are ellipses with fixed ellipticity 
€4 and position angle of the major axis 4>di & U centered on the photometric center (x c ,y c ). 
The intensity of the fcth elliptical isophote with semi-major axis is J^. 

The major-axes of all the ellipses make an angle 4>d to the x-axis of the image. We 
define elliptical coordinates (x e ,y e ) to be the position relative to axes centered on (x c ,y c ), 



3.1. 



Mathematical details 



We minimize the reduced xli which is defined in the usual way: 




(1) 
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and rotated to be aligned with the current estimate of major-axis position angle. Thus a 
pixel center at (x iy yi) has elliptical coordinates 

x e = (xi - x c ) cos (p d + (yi - y c ) sin <p d 
Ve = (Vi - Vc) cos (f) d - (xi - x c ) sin <p d . 

The semi- major axis, a iy of the ellipse passing through this point is given by 

2 



2 2 i 

a { =x e + 



Ve 



Since we adopt linear interpolation between the tabulated isophote values (the data do not 
merit anything more sophisticated), the two non-zero weights are 

Wk,i - 



W k +l,i = 



a k+l — a k 

cii — at 

a k+l — a k ' 



where {ak} are the semi-major axes of the K d ellipses at which the model disk intensities 
are tabulated. The intensity at the central point is I^ d and the model intensity for any pixel 
outside the the last ellipse is computed by linear extrapolation from the outermost pair. 

Our numerical task is to find the set of parameters {x c , y c , e d , <p d , Ik,d} that minimize Xr- 
For a given set of the parameters (x c , y c , e d , (f> d ), the values of {h,d} that minimize Xr can be 
obtained by solving the linear system 

K d / N \ N 

El \^ W j,i W k,i \ T \^ W j,i D i M\ 

k=i \i=i 1 / i=i 1 

obtained by setting the derivative dxl/dlj equal to zero. The search for the optimum choices 
for (x c ,y c , e d , (p d ) is by standard minimization; we adopt Powell's method (e.g. Press et al. 
1992). 

It is easy to include a bulge light component. We adopt a parametric Sersic bulge model 

l/n 



I s(r) = i s , exp < -b r , 



r , 

- 1 

r, 



(3) 



where n is the Sersic index, and I s $ is the intensity scale. The value of b n is set such 
that the effective radius, r e , contains half the light, and is defined by the implicit relation 
T(2n) = 2j(2n,b n ), with the T function and the incomplete gamma functions having their 
usual definitions. 
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We allow the bulge to be spheroidal, with the disk plane being the plane of symmetry. 
The parameter J s adds a single extra row to the matrix solution, whereas the bulge apparent 
ellipticity, e s , Sersic index, n, and effective radius, r e , add three non-linear parameters to the 
fit. We do not place bounds on these parameters, but sometimes find n increases without 
bound at every iteration, while J Sj o becomes very small. We found this behavior to be an 
excellent indicator of a bulgeless galaxy. 

Finally, we model the bar as an elliptical light distribution, with ellipticity and pro- 
jected position angle (fit, that can differ from those of the disk. Much of the foregoing can 
be applied with ease; we add an additional set of intensities {h,b\ to describe the bar, 
and all we need do is to extend the summation Ylk=i w k,dk such that K = Kd + Kb and 
{Ik} = {Ik,d} U {Ik,b}- We have used the same set of semi-major axes of the bar ellipses as 
for the disk, but our method could allow the disk and bar ellipses to be spaced differently. 

The Ik values are not guaranteed to be positive, and we regard significantly negative 
values as unacceptable, on the grounds that they would distort the relative light fractions 
in the different components. When the best fit model includes any significantly negative 
Ik values, we reset those Ik values to zero and re-fit. It is sometimes necessary to do this 
several times to ensure that the fitted model has no negative intensity values. (It is possible 
to enforce Ik > in each evaluation of xli but this results in discontinuities in the xl surface 
as ellipses are added and removed, which confuse the minimizing algorithm.) 

3.2. Procedure 

We selected a elliptical region of the image to be used in the fit. We first estimated the 
size and ellipticity of the galaxy by eye and then increased the semi-major axis by 30% and 
the axis ratio b/a by 10% to ensure that all pixels containing significant galaxy light were 
included. We masked out pixels containing light from foreground stars and other galaxies in 
order to eliminate them from the fit. 

We begin by fitting a disk only, and determine the four non-linear parameters x c , y c , €d, 
& (fid- Having found the best-fit center, we generally hold it fixed in all subsequent fits 
with additional components. Allowing the center to float in subsequent fits slows the fitting 
procedure by adding two extra parameters, but does not lead to a significantly improved xl- 
We then fit a bar and disk model, and a disk and bulge model, and finally a three-component 
model with a disk, bar and bulge. 

These steps are illustrated in Figure [1] for a sample case of the strongly barred galaxy 
g065 (PGC 035480). This Figure shows the fitted models, light profiles and residual maps 




Fig. 1. — An example to illustrate our procedure in a case with a clear bar, g065. The top-left panel shows 
the I-band image, with north up and east to the left, and the size in arcsec is marked by the scale bar. 
The other two panels in the top row show a disk-only model fitted to the data and the residuals when this 
disk is subtracted. Each pixel in the residual maps is weighted by its statistical uncertainty and the lighter 
shades indicate negative residuals and the darker positive, the elliptical outer edge bounds the fitted region. 
Masked pixels are unshaded. The next 3 rows show a disk+bulge fit, a disk+bar fit and our preferred three 
component disk+bulge+bar fit. The left-hand row shows the major-axis light profiles of the separate fitted 
components and the right-hand column the residual map in each case. Finally, the bottom row shows the 
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from each converged fit with the different set of components, together with the separate 
components of our final fit. 

The elliptical ring-like features in the centers of the residual maps (Fig. [U right hand 
column) for the two fits without a bulge component are artifacts of our linear interpolation 
scheme. Linear interpolation between our somewhat coarsely spaced ellipses does not allow 
the model to follow the steep bulge light profile all that well. However, this failing is of 
no consequence, since the steep light component in every galaxy is fitted by a Sersic bulge, 
leaving a shallower light profile to be fitted by the disk and bar, for which our interpolation 
scheme is entirely adequate. 

Ideally, we would like an objective test to determine whether the extra parameters of the 
additional components cause a significant reduction in the value of x% but we have not found 
anything satisfactory. We tried using the standard F-test, and also the Bayesian information 
criterion (Schwarz 1978; Mukherjee et al. 1998), but neither proved useful. The number of 
independent pixels in our images, after allowing for seeing, is perhaps one tenth of the actual 
number, (9(10 5 ), yet is still large compared with the number of parameters, (9(50), with the 
result that these formal tests nearly always indicate that the extra parameters are needed, 
even for a very small improvement to \ 2 T , such as might result from fitting part of a spiral 
pattern with the bar model. A further problem is that we do not attempt to model spiral 
arms, dust lanes, lop-sidedness, etc., and therefore we find typically 1.2 < Xr ~ 1-8, implying 
that all models are formally rejected with a high degree of significance, which undermines 
the rationale of these formal tests. We therefore rely on visual inspection of the residuals, 
after subtraction of the various models from the data, to determine the number of separate 
components in the galaxy. 

We have always employed equal spacing for the major axes at which the model intensities 
are tabulated, even though it is not a restriction of the procedure. Experience revealed that 
thirty ellipses to represent the disk was adequate; adding more did not reduce \l significantly 
for test galaxies, but Xr was significantly worse when many fewer were employed. We used the 
same semi-major axes for the bar, making a maximum of sixty ellipses when both components 
are fitted. The intensity profiles along a short bar are not determined very well by only a few 
ellipses, and with hindsight we might have spaced the bar ellipses more closely. However, 
the principal conclusions of this study are not affected by our undersampling of the bar light 
profile. 
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3.3. Discussion 

While the objective in our three component fits is to separate the disk from the bulge 
and bar, this did not always happen. In particular, the "bar" light fraction sometimes turned 
out to be an elliptical, presumably tri-axial, bulge misaligned with the disk major axis. In 
other cases with strong spiral arms and a weak or small bar, the fitting software sometimes 
identified these outer non-bar features as the principal non-axisymmetric component. These 
latter rare cases can be recognized by visual inspection of the model, and a new fit undertaken 
with fewer bar ellipses covering the inner disk only. 

The light profile of a "bar" that extends to the edge of the galaxy may not decline 
monotonically, but rise again well outside the visual bar where its major- axis crosses a spiral 
arm, for example. In order to find a more realistic bar, we repeat the fit with a smaller 
number of bar ellipses restricted to the inner region; we retain only those inside the first 
ellipse to fall below la of the noise in the first fit. The fits that result from this procedure 
seem reasonable in most cases, although the bar may be slightly larger than would have been 
estimated by eye. An alternative strategy would be to add a penalty to the xi function when 
the bar profile does not decrease monotically; we did not employ this strategy, however, since 
it would preclude the solution for the {Ik} by linear algebra. 

As our method relies on the shapes to separate the different components, there can 
be significant degeneracy between them if the shapes do not differ by much. For example, 
it becomes more difficult to separate bulge light from disk light when the galaxy is nearly 
face-on as both components are then assumed to have almost circular isophotes. In this 
circumstance, our method subtracts the best-fit Sersic bulge and assigns any remaining 
axisymmetric light to the disk. Such degeneracies become less problematic for more inclined 
disks. 

Here, we present our findings using the algorithm described above. We recognize that 
the light fractions in each component, and other parameters, may be systematically biased 
by our procedures. However, it is likely that every method to quantify bars suffers from its 
own idiosyncratic biases, and a method with a different set of biases is therefore useful. We 
are applying our technique to the OSUBSGS sample (Eskridge et al. 2002) that has been 
analysed by other methods, and will present comparisons elsewhere. 
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Fig. 2. — The first disk galaxy in our sampe to illustrate our fits. The full set of 68 disk fits is given in the 
on-line version. For each galaxy, the panels from left to right show: first, the cleaned I-band image on the 
left (North up and East to the left); second, the best fit model; third, the residuals in each pixel, weighted by 
its statistical uncertainty; and fourth, the surface brightness profile of each component in the fitted model. 
The scale is in I magnitudes per square arcsec for images taken in photometric conditions, and a simple log 
intensity scale otherwise. 



The apparent magnitudes of the galaxies observed during the first three photometric 
nights are listed in Table 1. The quoted values are isophotal magnitudes to the surface 
brightness 1=24 mag arcsec^ 2 . Comparison with existing measurements in the literature is 
possible in a few cases. Only two galaxies in our sample (g019 & g059) have isophotal I-band 
magnitudes given in the literature, and both agree with ours to within 0.05 mag. Since our 
B- & V-band exposures were deliberately rather short, our colors are not very precise. We 
found photographic magnitudes in the literature for 7 galaxies in the V-band and 9 galaxies 
in the B-band; our estimates agree with these published values to within their, rather large, 
errors. 



Of the 95 galaxies observed, we exclude six spiral galaxies (22, 28, 29, 57, 61 & 81) that 
had excessive light contamination from foreground stars, 2 galaxies with images that overlap 
another galaxy (64 & 67), and two galaxies with long tails indicating mergers or merger 
remnants (6 & 62). Throwing out these 10 cases reduces our sample to 85 galaxies. 



Our sample naturally contained elliptical galaxies. We identified and set aside the 
following 17 galaxies as ellipticals: gOOl, gOlO, g023, g027, g031, g038, g039, g044, g045, 
g060, g066, g069, g078, g089, g092, g093 & g095. 



4. 



Results 



4.1. Magnitudes 



4.2. Decompositions 
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The results of the fitting procedure are shown in Figure[2l From left to right, the I-band 
image, the best fitted model, residual map and light profiles of the fitted components are 
shown for all 68 disk galaxies. We give the fraction of total light in each component for 
the fit with the most components judged necessary to fit the image in Table 2. These light 
fractions report the total fluxes assigned to each component in our decompositions ignoring 
any flux in the masked pixels. Including masked pixels or extrapolating the bulge light to 
infinity would have hardly any effect on the light fractions. 

Uncertainties in the fluxes of the separate components were calculated by the Metropolis- 
Hastings algorithm (Monte Carlo Markov chain). We stopped the Metropolis evaluations 
once the changes in the 1 — a error bars dropped below 0.03% for an additional 20,000 
accepted points. We found that our parameters generally have small statistical uncertainties 
because of the large number of pixels. 

The first disk galaxy in our sample, g002 (NGC 3469), is a good example of a three- 
component model; the fitted disk, bulge (6.6%) and bar (28%), illustrated in the center 
panel, match the original image well. The residuals are dominated by the prominent spiral 
arms in the outer parts, but the regular pattern of residuals in the bar region indicates that 
the real bar is boxier than our fitted elliptical model. 

We note galaxy g043 (NGC 3479) has a smaller bar light fraction (6%) and is again 
fitted well, with the exception of the spirals. Galaxy g026 (NGC 4094) has the smallest 
identifiable bar light fraction (1.3%) of any in our sample, yet the bar is clearly present. A 
fit without it has a much higher Xr an d the bar stands out in the residuals. Galaxy g052 
(NGC 4924) is a good example of an unbarred galaxy that is fitted well by a disk and bulge 
only. 

4.2.1. Bulges 

Figure [3] shows the distributions of bulge light fractions in the barred and unbarred 
disks. Of the 68 disk systems, 12 have no bulge. The fraction of light in the bulge is 
generally moderate, and is greater than 20% in only 7 of the 68 cases, with the largest being 
31.5%. Apart from a slightly larger fraction of bulgeless disks in the unbarred sample, there 
does not appear to be a significant difference between the bulge light fractions in barred and 
unbarred galaxies in our sample. 

We find only a loose correlation, in the expected sense (e.g. Simien & de Vaucouleurs 
1986), between our estimated bulge light fractions and the Hubble T parameter. 
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Fig. 3. — The cumulative distributions of bulge light fractions in our sample of disk galaxies. The lower line 
shows the distribution in galaxies having identified bars, the upper line shows that in galaxies that appear 
to lack bars. 

4-2.2. Bar Fraction 

As is well known, identification of bars in disks becomes more difficult in more highly 
inclined galaxies. Among the 51 galaxies with apparent axis ratio b/a > 0.4, we identify 
bars in 24, i.e. almost half the cases. Although we have successfully fitted bars in a few 
more inclined disks, we do not find as large a bar fraction (4/17) in these cases - almost 
certainly due to this bias. Our 47% bar fraction in the less inclined galaxies is on the low 
side of recent estimates; the visual classification of the OSUBSGS sample by Eskridge et al. 
(2000) finds the bar fraction to be 70%, whereas the more quantitative criteria employed by 
Laurikainen, Salo & Buta (2004) and Marinova & Jogee (2006) yield bar fractions of 60% 
and 62% respectively in the same images. 

Our low bar fraction is all the more surprising since several of our barred galaxies have 
< 5% of the light in the bar; visual classifications may easily miss such small features. 
Furthermore, the bars in two other cases seem to be more triaxial bulges and may therefore 
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Fig. 4. — The fraction of bar light for 23 galaxies having RC3 types. The open circles are for galaxies 
classified SA, the triangles SAB, and the squares SB. Galaxies with no significant bar light are shown on the 
left with two separate values of the abscissae merely for clarity. The ordinate is the I-band magnitude where 
known, but is meaningless for the 4 galaxies below the horizontal dashed line, whose images were acquired 
in non-photometric conditions. 

have been classified visually as unbarred. Visual inspection of the images lacking bars by our 
criteria turned up just two possible cases that might have been visually classified as barred; 
in both cases, the possible bar is aligned with the inner end of a prominent spiral pattern in 
nearly face-on systems. However, it is quite out of the question that we could have missed 
as many as the 11 bars that would be required to bring our bar fraction up to 70% of the 
less inclined subsample of 51. 

It is possible that our sample of galaxies just happens to be deficient in bars; the 
discrepancy from Eskridge et al. is just 1.5y/~N, for our N — 51. Apart from small numbers, 
the discrepancy in bar fractions could be due to two possible systematic differences. First, 
Eskridge et al. used H-band images where the identifiable fraction of bars is higher than in 
visual bands; it is possible, but unlikely, that our I-band images also suffer from a similar, 
but lesser, bias. Second, the OSUBSGS is biased to higher surface brightness, since their 
sample selection criteria were B < 12 and D < 6'. Perhaps bars are less common in systems 
of lower surface brightness. 
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Fig. 5. — The cumulative distribution of bar light fraction in our sample of disk galaxies. 



4-2.3. Comparison with Visual Classifications 

Figure H] shows the Hubble classifications from the RC3 for 23 disk galaxies in our 
sample. We find six galaxies with SB classifications that have no bars, while we find bars 
with moderate light fractions for two galaxies that are classified SA. The original visual 
classifications were based on sky-survey images for the most part, which are of much lower 
quality than our data. 



4.3. Bar properties 

In all, 28 galaxies in our sample have an identifiable bar, for which the fit was signif- 
icantly improved by adding the bar component. As shown in Figure 0, the fraction of bar 
light has a broad distribution from 1.3% to 40% of the total light, with no preferred value. 

We estimate the bar length to be the semi-major axis at 10% of the maximum of the 
fitted bar intensity profile. Assuming the fitted disk is intrinsically round, we deproject the 
fitted bar length and ellipticity to find these quantities for a face-on disk. 

The bar light fraction correlates strongly with deprojected bar length, expressed as a 
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Fig. 6. — The deprojected bar length plotted as a function of the bar light fraction in our sample of disk 



fraction of -R25,/, as shown in Figure [6j Of the two points that lie below the trend, g013 
has a high inclination for which the deprojected the bar length is uncertain and g025 the 
"bar" may be a triaxial bulge. The single outlier above the distribution, g040, has a long 
bar in low surface brightness disk. As this correlation implies a somewhat constant surface 
brightness contrast for the bar, one worries whether the absence of large, low luminosity bars 
is simply a selection effect. After re-examination of the images, we think it unlikely that low 
contrast bars could have been missed, since the bars in all the cases on the upper edge of 
the distribution in Fig. [6] are still high-contrast features. 

Figure [7] shows the distribution of deprojected axis ratios, which has a broad range from 
0.35 < b/a < 0.8 (note that we plot 1 — b/a). Values of b/a closer to unity are excluded 
because such round bars would simply seem to be part of the disk. However, the lower limit 
is interesting - no bars in our sample are skinnier than about 3:1. Marinova & Jogee (2006) 
find a similar range of ellipticities from their different approach. It should be noted that this 
bound may depend on our use of simple ellipses; it seems unlikely that the distribution of 
values shown in this Fig. would be affected much by use of a more general bar shape, but 
the lower bound on b/a may be different. 

Figure [S] shows the disk light profiles, with the bar and bulge subtracted, of all the 
barred galaxies in our sample. The galaxies are ordered from top to bottom by increasing 



galaxies. 
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Fig. 7. — The distribution of deprojected bar ellipticities in our sample of disk galaxies. 

bar light frations. The light profiles of many disks decrease toward the center, which is 
mostly due to bulge subtraction. However, it is interesting that the bar component in a 
number of cases produces an additional inflexion in the inner disk profile, suggesting that 
the stars that make up the bar were taken from the radial range of the bar, and are not an 
additional feature "on top" of the disk. This feature is generally more easily recognizable in 
cases with large bar-light fractions towards the bottom of the figure. 



4.4. Colors 

We can estimate color differences between components by comparing component magni- 
tudes fitted to images taken in different filters. As our deeper I-band data is of higher quality, 
we adopt the best fit parameters for each component from the fit to the I-band image, which 
we then apply to the B and V images. We insist that the disk, bar and bulge have the 
same geometric parameters, and the same light profile in each of the different components. 
We therefore fit for the overall normalization of the light in each component in the B- and 
V-bands, using the projection parameters and light profile determined from the best fit to 
the I-band image. 

As shown in Figure M (upper), we generally find bars are slightly redder than the disk 
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Fig. 8. — The radial light profiles of the disks of all barred galaxies in our sample, scaled by the deprojected 
bar length. The different lines are arranged from top to bottom in order of increasing bar light fraction; they 
are shifted vertically by arbitrary amounts and shown with different line styles so that they can be easily 
traced. The vertical dashed line is drawn to show the bar length. The disks hosting bars that contain the 
larger light fractions (lower curves) generally show dips in their light profiles in the bar region. 

of the same galaxy and the bulges are significantly redder still, as shown in the lower panel. 
The small color differences between the disk and bar may indicate slightly different stellar 




B - I Disk B - I Disk 



Fig. 9. — Estimated colors of bars (above) and bulges (below) from our model fits. The line in each panel 
indicates equal color in each component. 

populations, or may simply be due to rather more dust in the bars. We have not corrected 
our plotted colors for internal extinction. 

Three galaxies (g008, g042, and g058) apparently have extremely red bulges. All three 
have rather faint bulges with < 5% of the total I-band light. One (g008) is quite edge-on, 
and all three appear to have considerable internal extinction that presumably contributes to 
the red colors. 



5. Noteworthy galaxies 

Two galaxies (g018 & g021) have prominent bars in otherwise low surface brightness 
disks. The fitted bar in g018 is 34% of the total light, although we feel this could be an 
overestimate. 

Our sample includes two apparently double barred galaxies: g019 (NGC 5878) and g070. 
The weaker bar feature in both cases, which is easily visible in the residual maps, appears 
close to perpendicular to the bright bar and seems to be of roughly equal length! The similar 
lengths of the bright and faint bars sets them aside from the "bars within bars" that are 
widely discussed (e.g. Erwin & Sparke 2002). 

In addition to g026, noted above, g004 and g098 are galaxies with extremely small bars. 
Our fitting procedure could detect these very small bars that may easily be overlooked in 
visual classification schemes. 

Four galaxies, g025, g030, g068 & g096, seem to have large intrinsic ellipticity differences 
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between the inner and outer disks. No one deprojection would yield a disk with near-circular 
isophotes; in all four cases the inner HSB disk appears to have a separate ellipticity from that 
of the outer LSB disk. The HSB inner disk in these cases does not correspond to a classical 
lens, since is has spiral arms in all four cases that are strongest in g030. The existence of 
spirals might suggest that the inner feature is the more nearly circular, which would imply 
a highly elliptic outer disk, or one with a very different inclination. 

Our algorithm selected a compromise ellipticity for the disk light of g025, which also 
has a small bar. Curiously, the position angle of the outer disk is within 20° of that of the 
bar. 

g030 was the most difficult galaxy in our sample to fit. Not only did the almost round 
inner HSB disk have strong spirals, but the more elongated outer disk was somewhat lop- 
sided. The fitting algorithm favored including the HSB inner disk as part of the slightly 
elliptical "bulge" and a more elliptical disk. There is some evidence of a bar in the visual 
image, but we were unable to obtain a 3-component fit that captured this feature - other 
decompositions were preferred. 

g074 and g094 are two moderately inclined disk galaxies (i ~ 45°) that both appear 
to have an excess of light on the major axis in two "blobs" symmetrically placed about the 
center. They do not appear as a classic bar, and may be a ring that is limb-brightened due 
to projection. 

6. Discussion 

We define disk size, -R/,25, to be the semi- major axis of 25th magnitude isophote in the 
I-band. Our definition of the bar length, as is the semi-major axis of the bar isophote at 
an intensity of 10% of its maximum. We find little correlation between Ri^ and in our 
models. This disagress with Erwin's (2005) finding that bar length scales with disk size, 
perhaps because we can find much smaller bars that may have been overlooked in previous 
work. 

We find no correlation with bar length and absolute magnitude, but find a good corre- 
lation with disk length and absolute magnitude as expected for systems of similar surface 
brightness. 

Ostriker and Peebles (1973) argued that the existence of a bar in a galaxy is determined 
by the dynamical importance of the dark halo; i. e. barred galaxies are predicted to have less 
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dark halo than have unbarred galaxieso Since it is well-established that dark matter varies 
inversely with galaxy luminosity, the Ostriker-Peebles stabilization mechanism indirectly 
predicts a decreasing bar fraction towards lower galaxy luminosity. We adopt a distance to 
every galaxy based on the SSRS redshift, using a Hubble constant of 70 km s _1 Mpc -1 . 

Since our sample is too small to look for a detailed trend of bar fraction with luminosity, 
we simply ranked our sample of nearly face-on galaxies by absolute luminosity. (Eight of 
the 51 barred galaxies in our sample were observed under non-photometric conditions and 
therefore lack I magnitudes.) The bar fraction in the fainter half (8/21) is slightly lower 
than that in the brighter half (13/22). However, a Kolmogorov-Smirnov test confirms that 
the difference in the distributions of absolute magnitudes between the barred and unbarred 
samples is of very low signficance. If our galaxy sample is representative, it would take a 
ten times larger sample for the minor difference between the luminosity functions of barred 
and unbarred galaxies to become significant. It is already interesting that our result is of 
low significance, since if bar-formation is inhibited by massive dark halos and if halo mass 
fraction varies strongly with luminosity, the luminosity functions of barred and unbarred 
galaxies should differ decisively. 



7. Conclusions 

We have described a method for separating the light of a disk galaxy into up to three 
separate components: a bulge, a disk, and a bar. The method also yields the radial light 
profiles of the disk and bar, as well as the axis ratio and length of the bar. The criterion for 
separating bar light from disk light is that the two ellipitical light components have different 
ellipticities and position angles that are radially invariant. We have not found an objective 
test to determine the number of components required by the data. 

We have applied this method to a sample of 68 spiral galaxies and shown that it works 
well in many cases, in that the separate components correspond to those that the eye identi- 
fies. In these cases, the light fractions and other properties of the separate components have 
moderate statistical uncertainties. 

The method is not completely objective, since it frequently fits features in galaxies that 
are not at all those intended. A triaxial bulge can be fitted as the bar, as can bright outer 



2 Athanassoula & Misiriotis (2002) reported an apparently contradictory result that bars can be excited 
by massive halos. However, such bars formed in their simulations have far too low a pattern speed to be 
consistent with bars in real galaxies (e.g. Aguerri et al. 2001). 
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spiral arms. Lop-sided or oval disks also generally cause problems. In most cases, the fit can 
be "nursed" to something close to the desired decomposition by restricting the radial range 
allowed for the bar fit, although in one case (g030) we were unable to obtain a satisfacory 
3-component fit despite repeated attempts with this strategy. 

Our procedure is still under development. One refinement that could be implemented 
quite readily is to allow for boxiness in the bar isophotes; the pattern of residuals after 
subtracting our model frequently indicates that a better fit would be obtainable by adding 
an extra bar shape parameter. Our method also has a tendency to taper the bar light out to 
larger radii than where the eye would judge the bar to end. It is fortunate that this weakness 
does not appear to compromise our estimates of bar length, since the low surface brightness 
of this excess bar light is below the 10% threshold that defines the bar light, but it may 
cause some overestimate of the bar light fraction. 

Although our statistical uncertainties are small, it is very hard to say much in general 
about systematic errors. These depend largely on how well our model fits the image - the 
fitted model matches the data well in some cases, but the light fractions in others are quite 
uncertain because the light distribution is more complicated than our simple 3-component 
model. Spiral arms are common, and several galaxies appear not to have intrinsically round 
disks; such features can, and clearly do, affect our estimates of the light fractions in the 
different components, but in ways that are impossible to quantify. 

This weakness aside, we find the fraction of barred galaxies in our sample is 47%, which is 
lower than other recent estimates. It is unlikely that visual classifications would overestimate 
the bar fraction, since our more objective scheme finds very small bars that could easily be 
overlooked by eye. Since galaxies were selected in an unbaised manner, its lower bar fraction 
could simply be due to our small sample size. 

The fraction of the total light that is in a bar ranges up to some 40%, with no preferred 
value. We have identified very small bars in a number of galaxies: the bar light in six 
cases is < 4% of the total, with the smallest identified fraction being 1.3%. The bar light 
fraction correlates with the relative size of the bar, suggesting a somewhat constant intensity 
contrast. 

The deprojected bar axis ratios range from 0.3 <b/a< 0.8. Intrinsically rounder bars 
would be hard to separate from the disk, but it is interesting that no bar was skinnier than 
3:1. 

We find some evidence that bars were formed from the disk, in that the more luminous 
bars weakly associated with disk light profiles that dip just interior to the bar end. We also 
find that bars appear to be redder than the host disk, but cannot say whether the color 
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difference is intrinsic or is due to extra dust in the bar. 

We find a slightly lower fraction of bars in the fainter half of our sample than in the 
brighter, but the difference of 5 cases out of 43 is barely significant. If disks are stabilized 
by halos, as originally proposed by Ostriker & Peebles (1973), then we should expect a lower 
bar fraction in low-L disks, which are known to be halo dominated. It is already interesting 
that the luminosity functions of barred and unbarred galaxies are so similar, but a larger 
sample of galaxies would provide a more decisive test of this hypothesis. 

We thank Roberto Abraham for providing some IDL scripts and both him and Mike 
Merrifield for helpful comments on a draft. We also thank an anonymous referee for an 
extensive and useful report. This work was supported by NSF grants AST-0098282 & AST- 
0507323 
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Table 1. Galaxy Photometry 
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1 Q 88 
10. oo 


i o fin 

1Z.0U 


1 Rn 

l.OU 


o 


.98 


oJOoO 


guis 


5591 


00907 


ppp okipka 

rblj U00004 


1 k a'? 

10. OO 


1 /I ^7 
14.0 ( 


i o /in 

1Z.4U 


o oc; 
z.zo 


1 


.28 




rrDI Q 

guiy 


5593 


00656 


lNJPP c;878 
IN KjK-j Oo ( O 


1 Q/l 

iz.y4 


1 1 70 
1 1 . 1 Z 


i n ^1 

1U.01 


1A 
Z.04 


1 


.32 




guzu 


5593 


00550 


ppp n^zL^87 

Jr VorO U04oo 1 


i p kp 

10.00 


i a on 

14. ZU 


1 A7 




Z.40 


1 


.35 




gUZl 


5497 


01342 


oa/t A qy Tino^c;oi i/mqic;i 

ZlVJTvoA. J 1UZOOOZ1-1440101 


1 0/1 
10. Z4 


1 /I 8fi 
14. oO 


1 Q IP 
10.00 


o qc; 

Z.00 


1 


.25 




guzz 


5506 


00659 


9A/TAQY T1 n^Ofi^Q 14^0^81 
ZlVJrioA. J lUOOZOOy-140ZOol 


1 K AP 
10.40 


1A 87 
14. ( 


1 3 8*3 
10. OO 


1 zLzL 
1.44 





.91 




gUZO 


5507 


01145 


PPP (\"XV77K. 
r uoo 1 1 


1 K A A 
10.44 


1 A C\P 
14. UO 


1 fil 
IZ.Ol 


^4 
Z.04 


1 


.23 




g024 


5514 


01149 


PGC 034225 


15.59 


14.24 


12.83 


2.33 


1 


.28 


L...P 


g025 


5514 


00991 


PGC 034459 


15.17 


14.19 


12.90 


2.04 


1 


.19 




g026 


5525 


00637 


NGC 4094 


13.08 


12.33 


10.94 


1.69 


1 


.05 


SXT6* 


g027 


5533 


01001 


2MASX J12295938-1437208 


15.96 


14.41 


12.78 


2.44 


1 


.27 




g028 


5533 


00831 




15.66 


14.59 


13.22 


2.44 


1 


.37 




g029 


5534 


01280 


PGC 042495 


15.55 


15.07 


13.91 


1.64 


1 


.16 


IB.9?/ 


g030 


5534 


01374 


PGC 042590 


16.19 


15.15 


13.93 


1.64 





.99 




g031 


5541 


00037 


PGC 043424 


13.87 


12.51 


10.93 


2.38 


1 


.34 


RLAR+P? 


g032 


5541 


00839 


IC 3831 


14.20 


12.90 


11.44 


2.30 


1 


.26 


PLXSO? 
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Table 1 — Continued 



Gal 


GSC 


Other Name 


B 


V 


I 


B-I 


V-I 


RC3 type 




5542 


00352 


NOG 4887 

IN VjrVj 4:00 1 


14. 


.26 


1 3 7fl 


19 74 


1 41 

X . T X 





.90 


T, A -l-P? 

Uii. |^ X . 








NQC, 4009 


12. 


.25 


1 1 ^9 

X X . d.Z 


m m 


1 Q1 

X . V X 


1 


.13 


OXJXVd 




5550 


01028 


TC! 4991 


14. 


00 


1 3 97 


19 14 

1Zi.1t: 


1 fi9 

X . \JZj 


1 


.00 


kJ-r\.XVdX ! 


pf)3fi 


5560 


00352 


NOG ^490 


14. 


.27 


1 3 


19 rn 


1 Q8 

X . C/O 


1 


.17 


n o* 

O . .d 


pf)37 


5576 


00288 


NOG ^664 


14. 


66 


1 3 73 

1U. I d 


19 41 

X — • JL X 


9 06 

Zj . WVJ 


1 


.22 




pf)38 


5593 


00419 


1NVJTVJ dOOU 


15. 


.54 


14 1 s ) 

It:, Xd 


19 ^f, 

XZj .dU 


9 4fi 


1 


.34 




guoy 


5593 


00376 


NOG ^88^ 

IN V.7 v_ OOOO 


14. 


.64 


1 ^ 94 

X d . Zjt: 


1 1 f,^ 

X X .Ud 


9 ^1 


1 


.36 




p-040 


5593 


00575 


9MASX 11^1^9060-1490317 


15. 


.39 


14 41 

It.t: X 


1 3 01 


9 04 

_ . v ' j: 


1. 


.21 




pf)41 


5505 


01314 


9MASX 110^31764-14^0^40 


16 


,42 


1 ^ 41 


14 19 

It . X <_/ 


9 07 


1. 


.14 




p-049 


5506 


01625 


9MASX T1 0^)^91 0-14^74^ 

1VX A kJ -/ V U X Wddd-Z X \J 14U 1 Tdd 


15. 


.46 


14 


1 ^ 4fi 


1 8Q 

l.Oi/ 


1 


.14 






5506 


01127 


NOG ^470 

IN V < V Ol I C/ 


14. 


.78 


1 ^ ^1 


1 9 1 Q 

IZj . X u 


1 Q 1 ! 
X . ZfO 


1 


.14 


SXR4 


p-044 


5507 


00560 


9MASX T1 lOO^I 8^-1 4^8901 

1VX kJ -/ Y O X X W <vd X Od 1t:uOaiW1 


15 


.84 


14 fifi 

14:.UU 


1 ^ 98 

lO .ZjO 


9 D8 

Zj. uo 


1 


.14 




p/045 


5507 


00627 


2MASX T11133590-1447377 


15. 


.75 


14 37 


12 87 


2.46 


1 


.30 




p-046 


6093 


00632 


9MASX 111^91816-1^00408 


15. 


.64 


14 79 

It:. 1 Zj 


1 3 43 

1U . t: - J 


9 07 


1 


.20 




p-047 


5525 


00508 


PGG 0^840fi 


14. 


.11 


1 ^ 91 


1 1 QQ 


1 8^ 

X . Od 


1 


.11 


PSRT^P* 

X OXJ X dX 


p-048 


5541 


00225 


ppp D44D91 


15. 


.44 


14 78 

It. 1 O 


1 ^ 7Q 

Id. 1 J 


1 4^ 

X .Hid 





.87 




p-040 


6111 


00332 


PGG n449fi7 

X V I V U 1 IZiU 1 


15. 


.17 


1 ^ 88 

X o.oo 


19 4fi 

IZj .tVJ 


9 0Q 


1 


.14 


SRT4* 

k J 1 J X t: 


ciWdW 


5542 


00818 


POO 04464 s ) 


15. 


.42 


14 39 


1 3 04 

Id . Wt: 


1 9 s ) 

X . u d 


1 


.18 


L -*/ 




5542 


00814 


PGG n44Q c )9 


14. 


.48 


1 ^ 7Q 

Xd . 1 c/ 


19 fi9 

IZj . 


1 f)Q 

X . dC/ 


1 


.00 






5542 


01198 


nhp 4094 

IN v_T V_ t:,^/ Zt: 


14. 


10 


1 ^ 9fl 


1 1 8 s ) 

X X .Od 


1 88 

X . oo 


1 


.13 


PSXSDP? 


p0^3 


5543 


01211 


9MASX T1 31 53736-1 459909 


15. 


.63 


14 53 


13 19 

Id . X U 


9 96 


1 


.26 










NGG ^07^ 

IN V_T V_ «JU I d 


13. 


.59 


1 9 fi4 


1 1 97 

1 l.Zi I 


9 1fi 

Zj . XU 


1. 


.31 


SRS 1 !?/ 

OXJkJd . / 


g055 


5550 


01415 


PGC 046523 


15. 


.47 


14.35 


12.95 


2.11 


1 


.18 




g056 


5551 


00322 


PGC 047994 


14. 


.82 


13.91 


12.34 


1.78 


1 


.05 




g057 






PGC 050081 


15. 


35 


14.21 


12.87 


2.47 


1 


.33 




g058 


5560 


00405 


PGC 050209 


15. 


.44 


14.43 


13.21 


1.66 


1 


.01 




g059 


5577 


00786 


NGC 5756 


13. 


.47 


12.56 


11.19 


2.01 


1 


.21 


PSBS4P/ 


g060 


5592 


01120 


PGC 053889 


15. 


.47 


14.14 


12.55 


2.53 


1 


.39 




g061 


6076 


00771 


IRAS 10495-1504 


15. 


.48 


14.86 


13.68 


1.57 





.99 




g062 


6077 


01161 


PGC 033374 


















g063 


6077 


01600 


2MASX J11050927-1521039 


15. 


.93 


15.25 


14.16 


1.52 





.92 
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Gal GSC Other Name B V I B-I V-I RC3 type 



g064 


6077 


00021 


PGC 034046 


15.48 


14.73 


13.69 


1.60 


0.95 




g065 


6085 


00590 


PGC 035480 












PSBR3 


g066 


6093 


00259 


PGC 037281 














g067 


6111 


00448 


PGC 043547 


15.68 


14.15 


12.40 


2.47 


1.36 




g068 


6111 


00192 


PGC 043625 


14.89 


13.93 


12.47 


2.07 


1.23 




g069 


6111 


00630 


NGC 4756 


13.85 


12.49 


11.04 


2.35 


1.29 


LXSO? 


g070 






NGC 4856 


11.90 


10.78 


9.42 


2.34 


1.31 


SBSO 


g071 


6112 


00735 


PGC 044701 


14.78 


13.92 


12.63 


1.80 


1.07 


IBS9?/ 


g072 


6112 


01482 


NGC 4877 


13.80 


12.52 


11.14 


2.16 


1.26 


SAS2* 


g073 


6113 


00188 


PGC 045958 












SXS8 


g074 


6113 


01709 


PGC 046436 














g075 


6122 


00596 


PGC 048102 


15.10 


14.58 


12.42 


2.25 


1.84 




g076 


6122 


00579 


PGC 048144 


16.26 


15.89 


13.34 


1.83 


1.65 




g077 


6123 


00184 


2MASX J13445714-1515151 














g078 


6153 


00233 


PGC 052002 














g079 


6154 


00667 


PGC 052840 














g080 


6168 


01004 


PGC 053634 












SBS7 


g081 


6062 


00839 


2MASX J10380158-1531057 














g088 


5525 


00508 


PGC 038496 


14.77 


13.72 


12.53 


1.80 


1.04 


PSBT3P > 


g089 






PGC 043646 


16.63 


15.25 


13.75 


2.45 


1.34 




g090 


6111 


01092 


PGC 043664 


15.10 


14.25 


12.84 


1.92 


1.24 




g091 


6111 


01258 


IC 0829 


14.62 


13.56 


12.35 


2.10 


1.12 




g092 


6111 


01164 


2MASX J12523303-1531010 


15.65 


14.38 


12.80 


2.30 


1.30 




g093 


6111 


00632 


PGC 043777 


15.88 


14.20 


12.83 


2.31 


1.09 




g094 






PGC 044471 


14.04 


13.14 


11.97 


2.00 


1.14 




g095 


6112 


00084 


2MASX J13071600-1543396 


15.80 


15.10 


12.84 


2.40 


1.83 




g096 






NGC 4984 


12.35 


11.26 


9.84 


2.22 


1.27 


RLXT+ 


g098 


6113 


00989 


PGC 046334 














gll2 


6095 


00798 


PGC 038529 


14.54 


13.30 


11.90 


2.63 


1.54 


SAS1? 


gll3 


6104 


01587 


2MASX J12390630-1610472 


15.38 


14.37 


13.09 


1.99 


1.13 




gl20 


6113 


01387 


PGC 046350 
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Table 1 — Continued 


Gal 


GSC 


Other Name B 


V I B-I V-I RC3type 


gl24 


6121 00362 


2MASX J13331239-1608216 




gl25 


6121 00343 


PGC 47717 


SXS2P* 



Note. - B,V,and I are R25 magnitudes. B-I and V-I are calculated at the 125 
isophote. 

b B uncertainty is +/- 0.048 mags 

V V uncertainty is +/- 0.035 mags 

: I uncertainty is +/- 0.046 mags 

*galaxy heavily contaminated by nearby star 



Table 2. Galaxy Decompositions 



Name disk % bar % bulge % e d <fi d (fib e« n r e 



g002 


65.31°;* 


oq 9+0.3 
^0-^-0.8 


6.618:1 


296+ 001 
u.zyu_ 012 


138.9!?;! 




U -" JO -0.01 


84.6!8i 


m+ - 01 

U - Ui -0.01 


0.83l8;8i 


3.418:1 


g003 


88.5l°; 5 




11.518^ 


n Q9+0.009 
u - oz -0.001 


124.618:? 






2^+ - 03 

u - zd -0.06 


0.5918:8! 


4.4i°:| 


g004 


74.5181 


I 7 i +0.4 
1 ' • i -0.4 


o cr+0.3 
0-3-0.4 


0-51818:88^ 


42.918J 


n q+0.01 

u -°-0.01 


7 7+1.9 
' • ' -2.1 


n7+ - 03 

u - u ' -0.03 


n 70+0.03 

u - ' °-0.05 


2 i +ai 


g005 


82.9l°i 


o 7+0.1 
°- ' -0.1 


i Q 4+ - 4 


n 901 +0.006 
u - zo± -0.008 


54.6l8j 


n c+0.02 

u - d -0.02 


no.9ll;5 


m+ - 01 

U - Ui -0.01 


1 m+ - 05 

i - uo -0.06 


0+0.2 

d - 8 -0.3 


g008 


97.9181 




2 l +ai 

^-O.l 


0.86718:881 


58.118;? 






0-4618:8^ 


o.63i8:8^ 


6-218:1 


g009 


QQ q+o-i 
yo.y_ 2 




6 l+°- 2 
°- 1 -0.1 


o.826i8:88? 


54.718;1 






0.2118:8^ 


0.9218:8^ 


5.2I8:? 


gOll 


73.5^q 6 




26.51?:? 


002+ 0009 

u - uuz -0.003 


91-71^ 






n 90+0.02 

u - zo -0.02 


1 iq+0-07 

i - iy -0.08 


1 2 9 +a4 
j-^-y-o.4 


g012 


77 7+O.7 
" •' -0.2 


16.31°^ 


6 +o.i 

°-0.5 


n 9C+O.OO2 
u - z " J -0.003 


148.518:1 


28 +a01 

u - zo -0.02 


39.71J 7 


m+ - 01 

u - Ui -o.oi 


o.85i8:8I 


3-918.! 


g013 


62.8iS:5 


qi cr+0.6 


5.6l? : ? 


o.73ii8:881 




0.2718:8^ 


102.61?;? 


21 +a02 

u - z,1 -0.06 


61 +a14 

u - ui -0.02 


2.618;? 


g014 


100 






n 1 70+O.OO5 
v - 1 1 °-0.007 


72.3l ; 2 












g015 


74.118:1 


19-6123 


6 4+ ai 


n 090+O.OO5 
u - ozo -0.006 


9 4+0.5 


o.56i8:81 


166.6181 


0-0218:8^ 


0.3118:8! 


5-718:1 


g016 


69.8l°; 2 




30.2I8:? 


0.78718:881 


ii8.7i8:l 






o.28i8;81 


1-6518:81 


9.418:1 


g017 


100 






n oQc+0.012 
u.zoo_ 006 


149.41?:^ 












g018 




341^ 


99+O.6 

zz -o.6 


001 +a002 
u.uui_ 001 


136.8l°:I 


n 90+0.02 

u - z,o -0.03 


42.4+0.9 


n 90+0.02 

u - zo -0.02 


o.75!8:8 2 


4.818:1 



Table 2 — Continued 



Name disk 



bar 



bulge 



4>d 



4>b 



n 



g019 


7318:1 


O+0.4 
O-O.l 


iq+ 01 

iy -0.4 


n fi4Q+ - 001 
u.u / ±y_ 001 


86.318:} 


fi+ - 01 

u -°-0.01 


101.418J 


o si +0 - 01 

u - oi -0.01 


1 91 +0.01 

i - zi -0.06 


9.6I8:? 


g020 


62.9+^ 


01 q+0.9 
oi - o -0.8 


5.818:1 


o.265i°:8? 7 


157l i 


^Q+ 01 
u - dy -o.oi 


4.618:1 


21 +a01 

u - zi -0.02 


fi7+ - 03 

u -° ' -0.04 


9 7+0.1 
z - ' -0.1 


g021 


60.61" 


20.31^ 


19.ll^ 7 


fl qqo+0.022 
U.OOO_q 012 


116.HL 


n c;9+0-06 

u - dz -0.02 


114 1 +L9 


n 99+O.O5 
u - zz -0.09 


n co+0.15 
u - DZ -0.06 


q c+0.9 

°- d -0.3 


g024 




40.3l?j 


n.4i8:^ 


n c;4+o.oo8 

u - d ^-0.006 


95.818J 


n 44+ - 01 
u -^-o.oi 


70.41?- 9 


1 Q+ 002 

u - ly -0.02 


n 7 cr+0.03 
u - ' d -0.03 


q+o-i 
o.y_ 2 


g025 


51.615:? 


30.71^ 


17.718:1 


n p-q+0.008 
u - do -0.008 


15.618:5 


26 +0 - 01 

u - zu -0.04 


183.411] 


0.0118:8? 


n 70+0-02 

u - ' y -0.04 


4.318:5 


g026 


94.2l8;i 




4.5i8: 5 7 


0.63918:885 


25.918:} 


n C.C.+0.02 
u - dd -0.02 


48.311J 


o.54i8:8^ 


n c;i + - 03 

u - di -0.04 


18.618:? 


g030 


88.41^ 




n.618:? 


o.i4l8:88^ 


5.9l 2 - 2 






o m+ 001 


n 77+O.O3 
U - ' ' -0.02 


418:1 


g032 


9018:1 




10I8I 


fl QQQ+0.002 


123I8:? 






n 90+0.02 

u - zo -0.01 


o.88i8:81 


4.818:? 


g033 


74.918:? 




25.118:1 


n 479+ - 005 

u -^' z -0.004 


113.218:1 






o.i4l8:8^ 


0.9418:81 


c+0.2 

8.O_ 02 


g034 


72 1 +ai 


i5.5l8:3 


12.418:1 


021 +0 ' 001 
u.uzi^q 012 


-9.il§:5 


o.65i8:8? 


23.1I8:? 


o.i8l8:8? 


1 01+0.01 

i -° i -0.09 


13-718.2 


g035 


97.61°:? 


i -°-o.i 


o.9l8i 


n ^99+0.003 
U.OZZ_ 003 


ioo.6i8:l 


o.68i8:8^ 


129l?j 


n7+ - 07 


0.6518:?? 


3 9 +0 


g036 


95.718:1 




4.418:2 


n fin4+ 001 
u.uu^^q 001 


133.318 






n 97+O.O6 
u,z ' -0.06 


1 9 c +0.08 
i - zd -0.07 


8.118:1 


g037 


81-6181 




18.418:1 


o c ;66 +0 - 002 


58.318:? 






0.0818:81 


1 14+0.03 

i - i ^-0.03 


5.718:? 


g040 


63.91};! 


1 o 7+O.6 
iz - ' -0.5 


23.31? 9 


u -^°-0.012 


74 7+0-9 

'^•'-0.8 


o.75l8:8? 


68.818:} 


0.1918:8? 


1 n3 +0 - 05 


4.318J 



Table 2 — Continued 



Name 


disk % 


bar % 


bulge % 








06 




n 




g041 


77 q+0.9 
1 1 .0_ 4 




22.t±8:S 


u.uoy_ .oi2 


74.818J 






o.i8i8:8^ 


0.92I8:? 4 


4-8181 


g042 


94.9l°;£ 




5.H5:S 


n c- 4 +0.002 
u -°^-0.003 


82.518:2 






o.i2l8:5 4 


n c+o.06 

u -°-0.04 


q 9+O.2 
°- z -0.1 


g043 


92.il8:S 


6 3 +a2 


1 6 +0 - 3 

i -°-0.3 


n 9 o 9 +o.oo5 
u.zoz_ 0004 


86.51™ 


41 + - 02 

u -^ i -0.02 


129.41H 


o.nl8:88 


0.5818:? 6 


q Q+0.5 

<J-O-0.6 


g046 


92.2l°j 




7 8+0.5 
'• 8 -0.6 


n qo+0.004 
u - oo -0.003 


196.518J 






o.o5!8: 06 


89 +0 - 07 

u - oy -0.08 


q c+0.2 
°-°-0.3 


g047 


57.4±g:i 


29.618:? 


1 '?+ - 1 
io -0.2 


n qcq+0.006 

U.OOO_g 


39.5l8i 


66+ 01 


44.9l8i 


o.o4i8:8^ 


OQ+0-02 
u - yo -0.02 


5.318:1 


g048 


69.718J 


30.318J 




n r 42 +0.009 
U.O^Z_ 012 


147.8l8i 


62 +0 - 01 

u - uz -0.01 


118.7181 








g049 


73.ii?:l 




4.2±8:i 


n qo+0.005 
u - oz -0.006 


"4.818:1 


u -^-0.02 


i5i.4lli 


o.H8:8^ 


n 79+O.O3 
u - ' z -0.05 


2.318:1 


g050 


100 






n 6O7+0-003 
u.uy i _o.oo3 


74.318^ 












g051 


100 






U. z ±ZZ„ 005 


101.8181 












g052 


95.218:1 




4.818:1 


f) 1Q1 +0-003 

u.±y±_ 003 


32.718J 






o.o8l8:8? 


o.62i8:81 


3.618:1 


g053 


Q9 q+0.2 




7 i +0.2 
' • i -0.2 


0.63918:88^ 


77.5l8i 






2 c i +0 ' 06 


n 77+O.O6 
u - ' ' -0.03 


c c;+0.2 

O.O_ 2 


g054 


92.7;°i 




7.318:1 


0.84418:881 


120.318:1 






47+ - 01 

u -^' -0.01 


o.88i8:8? 


9.418:1 


g055 


OCT O+0.4 

oO.O_ 4 




14.7181 


n 4Q7+0-0 07 


209.918:3 






09 +0 04 

u - uy -0.03 


i 91 +0.06 
i - zl -0.06 


4.2181 


g056 


64.9±g : | 


23 7 +1 - 4 

z "->- ' -0.9 


ii-518:2 


n nsi +0-016 

U.U81_o.oo6 


-38lf 


42+ - 02 

u -^ z -0.02 


121.3181 


n 99+O.O4 
u "^-0.01 


o p;s+ - 06 

u -°°-0.02 


9 7+0.1 
z - ' -0.1 



Table 2 — Continued 



Name 


disk % 


bar % 


bulge % 




4>d 




fa 




n 




g058 


97.31°;! 




O 7+0.1 
z - ' -0.1 


n 1 4 +o.oo9 

u - i ^-0.012 


57.9l?;J 






m+ - 01 

u - Ui -o.oi 


u - ud -0.06 


3.918:1 


g059 


92.il8:i 




7 q+o.1 
' - y -o.i 


f) coi +0.001 

u.ooi_ 0001 


48.5181 






1 3 +0 - 01 

u - io -0.01 


1 1 c+0.01 
i - id -0.02 


4.4181 


g063 


100 






n RQo+o-ooe 
u.uoy_ 008 


141181 












g065 


79.2±° : | 


13-3185 


7 c+0.2 

' -°-o.i 


n 241+ - 007 

u - z ^ ± -0 .005 


13.115:1 


0.5918:8^ 


69.5l8i 


ni+ - 01 

u - u± -o.oi 


o.68i8:8? 


4.9181 


g068 


96.7l°:l 


3.318:1 




266+ ' 004 
u.zuu_ 004 


192.9181 


0.4518:81 


16.51H 








g070 


49.318:1 


19.21°:? 


31.5181 


u -°+' -0.001 


51.4181 


o.48i8:81 


qpr+0.3 
od -0.2 


n 0+0.01 
u -°-o.oi 


1 qq+0.01 

i -°°-O.01 


17.2I8I 


g071 


100 






n fi 77+o.oo3 

u - u ' ' -0.003 


i4o.9i8i 












g072 


yd -0.1 




c+0.1 
°-0.1 


n c;79+0.001 
u - d ' z -0.001 


82.7181 






1 7 +0 - 02 
V - L 1 -0.01 


0.9218:81 


°- y -o.i 


g073 


99.318:1 




0.718:1 


001 + 001 
u.uui_ 001 


80.11& 






1 c i +0 - 08 

u - id -0.06 


o.67i8:8? 


4.718:^ 


g074 


90.718:1 




9.318:1 


0.71518:881 


165.7181 






u - od -0.01 


o.sgi^ 


6.318:1 


g075 


92 




7 c+0.2 
' - d -0.2 


o.448i8:88^ 


i5o.8l8i 






0.0418:8^ 


o.88l8:8^ 


4.H8:1 


g076 


95.9±?:i 




4.H?:I 


n i r 1 +0.012 
u - idi -0.015 


ioo.4il:| 






0.1818:8! 


0.5918H 


3-9111 


g077 


100 






n -| 4 r+0.015 
U.14O_ 0024 


83.31^ 












g079 


77.2181 


15.618:5 


7 q+0.8 
' -°-0.8 


n O97+0.009 
U.OZ, i _ .006 


154.4181 


0.4318:8^ 


180.6l8i 


0.1718:8^ 


0.8218H 


q 7+0.2 
°- ' -0.4 



Table 2 — Continued 



Name disk % bar % bulge % 



4>d 



4>b 



n 



g080 


100 






U.OOO_o.ooi 


143.8±5. 6 












g088 


98.2±g;| 


1 8 +0 - 2 

L -°-0.3 




o.o7H8:8f 


39.4181 


0.6218:8^ 


27.iil:l 








g090 


9i-2±8:i 




o O+0.9 
O-O-0.6 


o.369l8:88^ 


109.818:1 






u - zi -0.08 


o.74i8:8! 


a 9+O.3 


g091 


64.91J:! 


14.21°;? 


20.91^ 


n iqq+ 01 


156.41^ 


0.2918:8^ 


ii6.9l8:| 


24+ - 01 
u - Z4 -o.oi 


o.6H8:8^ 


418:1 


g094 


79.il8:i 




20.918:1 


fiqi+ 001 


178.518:1 






n o+o.oi 
u -°-o.oi 


0.841811 


4.518:1 


g096 


98.3181 




1 7 +ai 
x - ' -0.1 


n i qa+o.ooi 
u.io^_ 001 


oq O+0.2 






ni+ - 01 


0.3918:81 


2.718:1 


g098 


84.418;! 


3.918J 


1 1 7 +0 - 2 

ll -< -0.8 


U.4DO_ 003 


ii5.6i8i 


0.4718:81 


147.71H 


1 7 +0 - 04 

V - L ' -0.03 


0.911!:;;:; 


4.518:1 


gll2 


84.8181 




i c O+0.6 
±o - z -0.2 


0.50318:88^ 


106.3181 






n i 0+0.02 
U.18_ 002 


0.8518:81 


4.918:1 


gll3 


89 8+ 8 - 1 
°^-°-0.2 




i7.2l8i 


0.08218:81 


180.515:1 






0.031818 


0.67181? 


a 9+15.6 
4 - z -13.2 


gl20 


98.518:1 




i-5l8:J 


0.58618:881 


62.818 






0.34181^ 


o.58i8:8s 


3.4l? 3 4 2 


gl24 


100 








150.4l|8 












gl25 


94.il8:i 




5.918:1 


o c ;o6 +0 - 001 

U.OUO_ .oo6 


153.1181 






o.i9l8:81 


Q8 +a03 
u - yo -0.05 


4.5181 



a Disk, bar, and bulge refer to the fraction of light in each component. 



